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ABSTRACT 


A  computer  program  is  presented  which  solves  the  search 
planning  problem  for  survivors  at  sea.  The  program  is 
designed  to  utilize  weather  data  as  compiled  by  the  United 
States  Navy  at  its  Fleet  Numerical  Weather  Central/ 
Monterey/  California. 

An  investigation  is  also  made  into  the  search  criteria 
used  by  the  United  States  Coast  Guard  in  its  planning 
procedures.  Guidelines  are  given  for  the  use  of  the  square 
search  and  the  Sector  search. 
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I.  INTRODUCTION 

A.  THE  COAST  GUARD  AND  SURVIVOR  SEARCH 

Survivor  search  Is  defined  as  the  process  of  planning 
and  executing  a  search  for  survivors  at  sea,  and  the  Coast 
Guard  usually  accomplishes  this  In  the  following  steps: 

1.  Initial  Notification 

This  stage  Involves  gathering  all  available  data 
about  persons  or  craft  In  distress.  Including 

characteristics  of  the  survivor  craft,  number  of  survivors, 
etc. 

2.  Initial  Distress  Position 

In  this  stage,  the  search  planner  puts  together  all  the 
available  data  concerning  the  survivor  and  makes  a  judgement 
as  to  what  the  most  likely  position  of  the  survivor  was  at 
the  time  the  craft  seased  operating  normally. 

3.  Drift  .Calculations 

The  survivor's  suspected  motion  plotted  assuming 
drift  In  the  open  ocean,  and  a  position  is  calculated  for 
some  future  time,  usually  the  planned  beginning  of  search 
efforts. 

4.  Calculation  of  the  Area  to  be  Searched 

An  estimate  is  made  of  the  accuracy  of  the 
survivor's  position  as  calculated  in  step  3  and  the  accuracy 
of  navigation  of  the  search  craft.  An  area,  called  the 
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Search  Area#  Is  defined  based  on  these  errors 


5.  Allocation  of  Search  Units 

A  determination  Is  made  of  the  resources  (ships, 
planes,  boats,  etc.)  to  be  used  In  the  search,  and  the  way 
In  which  these  resources  will  be  used. 

6.  Execution  of  the  Search 

Searching  Is  accomplished  according  to  several 
established  search  patterns,  most  of  which  are  different 
arrangements  of  a  series  of  parallel  paths  within  the  Search 
Area. 

7.  Termination  of  the  Search 

Either  the  survivors  are  located  or,  after  a  certain 
amount  of  searching  Is  completed  without  success,  the  search 
Is  terminated.  This  stage  Includes  reports,  and  the 
evaluation  of  search  effort. 

B.  DOCTRINE 

The  guidance  for  all  these  activities  Is  contained  In 
Coast  Guard  publication  CG-308,  the  National  Search  and 
Rescue  Manual,  (Ref.  1.).  Chapters  3  through  8  of  this 
manual  contain  the  Instructions  for  planning  and  conducting 
a  search  for  survivors  at  sea. 
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C.  PURPOSE 

This  thesis  accomplishes  two  things: 

1.  Provides  a  computerized  solution  to  step  3  of  the 
search  planning  process#  calculating  survivor  craft  drift# 
and 

2.  Investigates  the  criteria  used  for  allocation  of 
resources  to  a  search#  steps  4#  5#  and  6. 

It  is  Intended  that  the  results  be  understood  by  the 
search  planner#  so  that  they  may  be  of  some  use  to  him. 
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A.  OBJECTIVE 

The  major  objective  in  designing  a  computer  program  for 
search  planning  was  to  have  a  fast,  accurate  means  of 
obtaining  an  estimate  of  a  survivor's  oceanic  position  at 
the  time  of  the  beginning  of  a  search. 

B.  THE  MANUAL  METHOD 

in  order  to  understand  the  program.  It  Is  first 

necessary  to  understand  how  drift  plotting  Is  done  manually. 

i.  Hamate 

A  fishing  vessel  Is  believed  to  have  sunk  at 
latitude  40-42. ON,  and  longitude  52-13. 0W,  at  0900z  on  22 
June.  The  probable  error  of  this  position  Is  estimated  to 
be  30  miles.  Although  the  standard  measure  of  variability  Is 
the  standard  deviation,  all  error  estimates  made  In  the 
search  planning  process  are  In  terms  of  probable  error.  It 
will  be  demonstrated  that  this  error  Is  defined  by  the 
statistical  term,  "probable  error".  Survivors  are  believed 
to  be  In  rubber  rafts. 

The  proper  Coast  Guard  authorities  are  notified  at 


2100Z  on  22  June/  and  a  search  plane  has  been  dispatched.  It 
is  due  to  arrive  on  the  scene  at  0100Z  23  June.  The 
aircraft's  navigation  error  (atain  a  probable  error)  is 
estimated  to  be  10  miles. 

2.  (lathering  Data 

The  first  task  of  the  search  planner  Is  to  estimate 
wind/  wind  current*  and  sea  current  at  the  position  of  the 
sinking.  Mind  Is  obtained  from  observations*  weather  maps* 
Navy  Fleet  weather  services*  or  the  U.S.  Weather  Bureau* 
whichever  Is  available.  Wind  current  is  computed  from  the 
wind  by  means  of  a  graph  given  In  the  National  Search  and 
Roscue  Manual.  Average  Sea  Current  is  obtained  from  ocean 
current  chr.rts  or  pilot  charts.  These  sources  for  average 
sea  current  are  compiled  from  observations  made  over  several 
years.  (Thus  this  data  may  not  be  accurate  for  the  time  and 
place  of  concern*  but  so  far  they  have  been  the  only  data 
available  to  the  search  planner.) 

3.  Drift  Pic ‘-tine 
e.  Leeway 

Leeway  Is  the  motion  of  the  si.rvlvor  craft  caused  by 
wind.  For  the  examole  rase  of  a  life; aft*  two  graphs  are 
provided  In  the  SAR  Manuel  for  computing  ’eeway.  One  graph 
Is  for  when  e  drogue  Is  attached  to  the  raft  (a  drogue  Is  a 
type  of  sea  anchor  designed  to  retard  llferatr  leeway)  and 
the  other  Is  for  when  e  drogue  Is  .ot  attached,  if  it  is  not 
known  whether  e  drogue  Is  attaches,  both  conditions  are 
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considered  separately. 

b.  Total  Drift 

Once  leeway#  wind  current#  and  average  sea  scurrent 
have  been  determined#  they  are  added  vectorial ly  to  produce 
total  1  I feraf t  drift. 


This  calculation  can  be  done  for  Intervals  of  any 
length  -  one  hour#  three  hours#  twelve  hours  -  but  the 
Interval  most  commonly  used  Is  the  interval  between  wind 
observations.  That  Is  If  wind  is  available  for  0900Z,  1200Z# 
1500Z#  etc.#  the  search  planner  would  probably  compute  total 
drift  In  three  hour  Intervals#  from  0900Z  on  the  22nd  to 
0100Z  on  the  23rd.  Each  time  total  drift  Is  calculated 
(another  probable  error)#  12.5%  of  that  distance  Is  retained 
as  "drift  error"#  and  the  drift  errors  are  summed  to  obtain 
total  drift  error. 

Example:  In  this  manner#  the  position  of  the  example 
llferaft  Is  calculated  to  be  40-58. ON#  51-55. 0W#  and  the 
drift  error  Is  4.0  miles. 

c.  Minimax  Plotting 

If  it  Is  unknown  whether  a  drogue  was  in  use# 
two  plots  are  made#  and  the  final  ("datum")  position  is 
defined  to  be  midway  between  the  two  positions  calculated. 
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The  drift  error  is  a  function  of  the  two  calculated  drift 
errors  and  the  distance  between  the  two  positions.  This 
process  of  plotting  two  positions  Is  called  Mlnimax  plotting 
in  the  SAR  Manual. 

d.  Error  of  Posit  lor 

Once  the  datum  has  been  computed.  It  Is  then 
necessary  to  find  total  probable  error  of  position,  c,  where 

de  ■  total  drift  error 

X  ■  error  in  initial  position  of  survivor  craft 

Y  ■  nav  1  gat  1  on  error  of  search  craft 

c  =  \[7rrTxrZ  y * 

For  the  example,  c«  ^16+100^900  *^1016  •  31.9  miles. 

e.  Safety  Factor 

To  determine  the  size  of  the  area  to  be 
searched,  a  search  "radius",  R  Is  found  by  multiplying  c  by 
a  safety  factor,  depending  on  which  search  in  a  sequence  is 
being  planned.  For  searches  one  through  five,  the  respective 
safety  factors  are  1.1,  1.6,  2.0,  2.3,  and  2.S.  The  example 

of  the  sunken  fishing  boat  is  a  first  search,  so  R  ■  1.1  x 
lk. 7  -  16.2  miles. 

The  area  to  be  searched,  then.  Is  a  square  with 
sides  ■  2xR,  centered  on  datum: 
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f.  Sweep  Width 


The  standard  measure  of  the  effective  search 
width  of  a  search  craft  is  called  "sweep  width".  The  value 
of  the  sweep  width  Is  obtained  from  a  table  In  the  Search 
and  Rescue  Manual ,  based  on  the  type  of  target,  search 
altitude,  and  meteorologl cul  visibility.  The  sweep  width 
thus  obtained  is  corrected  for  the  effect  of  st-a  state  by 
multiplying  by  a  "whltecap  correction  factor"  which  is 
tabulated  for  various  wind  speeds. 

C.  THE  COMPUTER  PROGRAM 

1.  QhjlfcLlM 

The  computer  program  accomplishes  the  objective  of 
providing  an  estimate  of  the  survivor's  position  by 
automating  the  steps  of  the  manual  method  of  drift  plotting. 

*.  Sflurci  af  Miiihir  Bali 

In  choosing  e  source  of  weather  data  for  the  computer 
program,  three  criteria  were  considered:  a.  aval  lability,  b. 
computer  access  tin',  and  c.  accuracy.  Th«  final  selection 
for  data  source  was  the  Navy  fleet  Numerical  Heather 
Central,  Monterey. 

a.  Availability 

FNWC  Heather  data  is  compiled  twice  dally,  at 
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OOOOZ  and  1200Z,  and  is  directly  available,  either  by  data 
link  or,  for  testing  purposes,  on  magnetic  tape.  For  use 
with  this  program,  the  data  is  stacked  chronologically  or  a 
master  tape. 

b.  Computer  Access  Time 

With  the  exception  of  the  mounting  of  the  master 
weather  tape,  all  weather  data  handling  is  eliminated. 
Access  speed,  then,  is  limited  only  by  the  speed  of  the  tape 
transport . 

c.  Accuracy 

Although  the  data  used  is  intended  for  wide  area 
coverage,  its  probable  error  is  less  than  the  probable  error 
assumed  when  computing  drift  error  (12. 5%).  Although  there 
will  be  occasions  when  local  weather  reports  are  available 
and  may  be  more  accurate,  the  bulk  of  search  planning  is 
done  under  conditions  of  incomplete  weather  analysis  and 
forecasts,  and  when  gross  estimates  must  be  made.  Since  the 
weather  analyses  for  a  particular  time  period  and 
geographical  location,  once  assembled  for  this  program,  are 
fixed  parameters,  repeated  computations  of  the  same  search 
problem  will  always  produce  Identical  answers.  This  permits 
a  rigorous  evaluation  of  the  data,  methods,  arid  assumptions 
used  in  search  planning  by  a  systematic  testing  procedure. 

3.  DaU.  Description 

The  weather  data  are  produced  by  FNWC  for  3,969 
geographical  points  in  the  northern  hemisphere.  These  data 
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points  are  intersections  of  a  63  X  63  square  grid 
superimposed  on  a  polar  stereographic  map  projection.  Grid 
point  <32, 52 )  is  at  the  north  pole,  and  the  vertical  line 
< i “32)  is  coincidental  with  the  10  W  and  100  E  meridians.  A 
sample  projection  is  shown  as  Figure  1. 

The  polar  stereographic  projection  provides  a 
plotting  reference  which  shows  an  entire  hemisphere  with  the 
least  possible  distortion  (Ref.  6),  and  the  scale  is  true  at 
60  degrees  North  Latitude.  Conversion  from  latitude  (L)  and 
longitude  (*)  to  grid  coordinates  is  accomplished  by  the 
following  formulas: 


<’  -  32  +  31  20‘J 


<  -  SinL  * 
*  -  *‘nL 


Ct  b  ( *  *  ) 


(1) 


J  W  -  Jl  2c*f 


*  -  k»  r»  (_  1  ‘  > 

- —  3  in  (X  ♦  »<  ) 

/  .  Sin 


(7) 


For  the  esample,  I  •  31.001,  and  j  •  19.0163. 

The  date  Is  stored  In  "fields"  of  3,919  words 
(numbers),  each  word  being  some  component  of  atmospheric  or 
oceanic  conditions  at  Its  corresponding  (i,j)  grid  position. 
The  fields  required  for  direct  Input  to  this  program  are: 
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U  Curr 


U-  component  of  surface  current 


V  Curr  V-  component  of  surface  current 

U  Narine  U-  component  of  surface  wind 

V  Marine  V-  component  of  surface  wind 

where  the  U  direction  is  parallel  to  the  "I"  or 

horizontal  ails  and  the  V  direction  Is  parallel  to  the  "J" 
or  vertical  ails,  in  the  FNWC  grid  system. 

The  Current  information  as  provided  by  FNWC  is  total 
surface  current*  including  wind  current. 

Fleet  Numerical  did  not  regularly  produce  u  Marine 
and  V  Marine  wind  fields*  so  these  fields  are  computed  by  a 
separate  computer  program  provided  by  the  Weather  Facility 
using  the  following  data  fields  as  Input: 

0  1000  Height  of  the  1000-mto  pressure  line 

T  Air  Air  temperature  at  the  surface 

T  Sea  Sea  temperature  at  the  surface 

0.  CALCULATING  OR  I  FT  WITH  THE  COMPUTER 

1  •  Leeway 

In  addition  to  the  llferaft  drift  tables  given  in 
the  Search  and  Rescue  manual,  there  are  other  ways  to 
predict  the  leeway  of  a  survivor  craft.  Most  of  these 
methods  Involve  using  a  percent  of  the  wind  speed  to 
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determine  the  leeway  speed,  *nd  some  deflection  in  degrees 
to  determine  the  direction  of  leeway  off  the  wind  line.  For 
example,  a  boat  type  may  be  known  to  drift  in  a  direction  40 
degrees  off  the  wind  line,  and  at  a  speed  equal  to  4$  of  the 
wind.  To  accomodate  these  methods  as  well  as  the  llferaft 
drift  table  from  the  SAR  Manual,  several  options  are 
available  to  the  user  of  this  computer  program. 

2.  Options  for  Leeway 

a.  A  function  approximating  the  ltferaft  drift 
tables  in  the  National  SAR  Manual  has  been  provided; 
drifting  speeds  for  wi th  drogue  and  wi thout  drogue  are  both 
aval lable. 

b.  Leeway  speed  can  also  be  specified  as  a  percent 
of  wind.  In  whole  numbers,  from  2\  to  9i. 

c.  For  each  of  the  leeway  velocities  in  b.,  a 
deflection  off  the  wind  line  of  plus  or  minus  some  angle,  in 
tens  of  degrees,  from  0  to  90  degrees,  may  be  specified. 


d.  A 

mi ni mum-max  1  mum 

1  eeway 

speed  option  wi th 

the 

minimum  being 

1%  of  the  wind. 

and 

the  maximum 

bel  ng 

some 

whole  number 

percent  from  2% 

to 

9%  of  wind 

speed. 

i  s 

aval lable. 

e.  No  leeway  at  all  (as  in  the  case  of  a  man  in 
the  water  or  some  other  almost  completely  submerged  object) 
may  be  specified. 

3.  Adjustment  of  Input  Data 
a.  Dates 
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All  dates  are  converted  from  month  and  day  to 
Julian  dates,  to  permit  reference  to  weather  data  and 
adjustment  of  the  program  "clock". 

b.  Latitude  and  Longitude 

All  Initial  positions  are  converted  to  FNWC 
(l,J)  coordinates. 

c.  Leeway  Options 

Leeway  options  are  entered  by  means  of  j 
two-digit  code,  which  Is  converted  to  speed  and  direction  of 
drift  within  the  pregram. 

d.  Weather 

Wind  and  current  are  provided  in  units  of 
cm./sec.  and  are  scaled  In  the  program  to  nautical  miles  per 
hour.  They  are  further  scaled  to  decimal  fractions  of  the 
FNWC  grid  distance. 

k.  Interpolation  of  Weather  Data 

a.  Interpolation  for  Position 

For  each  position  In  (l,j)  coordinates  there  are 
four  neighboring  grid  Intersect  Ions  for  which  weather  Is 
tabulated,  (sample: 

CHID  POSITION:  l*59. 007IT 

J*19. 01*15 

0*10  POINTS:  .  (51,10)  .  (59,10) 

.  (51,19)  .  (59,19) 

All  weather  data  are  interpolated  linearly  among 
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these  four  neighboring  points.  Fir  example,  if  (J  Current  at 
0000Z  at  grid  points  (l,j)  is: 

( 31, 20)  -  SO.O  (  59,  20)  -  *9.0 

(31,19)  ■  *2.0  (59,19)  •  *0.0 

then  for  l»3l,  j"19. 01*13, 

U  Current  •  *2.0  ♦  . 0 1  %  1 3  x  (S0-*2)  -  *7.117 
and  for  I "39,  a  similar  computation  results  In  U 


Current  •  *0.12717. 

U  Current  for  the  position  is  then, 

*2.112  ♦  .00717  «  (-1.9IS)  •  *7.097  cn./sec. 


b.  I  nterpolat ion  for  Time 

For  each  hour  in  the  computation,  there  are  two 
observation  times  which  are  used  to  supply  weather  data: 
00002  and  1200Z,  one  of  which  will  be  prior  to  the  hour,  and 
one  which  will  be  after  the  hour.  A  linear  I nterpolat Ion  for 
oosltlon  is  done  for  each  observation,  and  then  these  two 
values  are  interpolated  linearly  for  time.  Tor  example,  if 
the  U  Curr  for  17002  is  calculated  as  above  and  found  to  be 
37.25  cm. /sec.,  the  Interpolat Ion  for  time  is  accompli 
as  follows: 


*1.09  ♦  9/17  x  (-1.1*)  •  59. *f  cm. /sec. 
c.  The  Cloc* 

A  complete  set  of  I  nterpolat  Ions  is  done  for 
each  position  and  for  each  weather  parameter  at  the  start  of 
the  drift  calcu I  at  I ons,  and  each  position  is  advanced  for  an 
hour's  worth  of  drift.  This 
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interval  is  used  as  a 


compromise  between  accuracy  requirements  and  computer  time. 
If  an  Interval  much  smaller  than  one  hour  were  used,  a  large 
Increase  in  computer  time  would  be  required.  A  problem 
which  would  take  20  seconds  at  one  hour  Intervals  would  cost 
almost  20  minutes  of  computer  time  if  one  minute  Intervals 


After  the  one  hour's  worth  of  drift  is  computed, 
and  the  positions  updated,  the  clock  Is  advanced  one  hour, 
and  the  process  Is  repeated  over  and  over  until  the  end  of 
the  problem. 

s.  Drift  Platting 

Drift  calculations  are  applied  to  the  (i,J) 
coordinates  of  the  survivor  craft  plot.  That  Is,  the 
U-component  of  drift  Is  added  to  the  "I"  coordinate  of  the 
survivor  position  and  tha  V -component  Is  added  to  the  "J" 
coordinate. 


examples  The  Interpolated  values  of  current  and 
wind  for  OtOOZ  on  22  June  are 


Pan 

■meter 

cm. /sec. 

ml les/hour 

U 

Curr 

♦  Sf.klO 

♦  .711 

V 

Curr 

♦  51.121 

♦  .702 

U 

Mar  1 ne 

♦250.71 

♦  k.17 

V 

Marine 

•ISO. 25 

-11.51 

tha 

wind 

speed  then  equals  w4.67**ie.sy*' 

■17.25 

nautical  miles/hour 

•  For  a  llfaraft 

without  drogue. 

this 

wind  spaed 

produces 

a  leeway  of  0. 

9k  miles  per  hour  as 
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f o 1 1 ows : 

Leeway  «  .  226  x  ^ 17. 25  »  0.94 

This  formula  is  an  approximation  to  the  iiferaft 
drift  table  given  in  the  Search  and  Rescue  Manual. 

Since  arctan  (4.87/-16.55)  equals  163.6  degrees, 
this  angle  is  a  measure  of  the  wind  angle  (measured 
clockwise  from  the  vertical  (In  the  FNWC  system): 

V 

If  leeway  deflection  is  specified  by  the  user,  it  is 
applied  to  this  wind  angle,  before  the  U  and  V  components  of 
leeway  are  calculated.  The  reduction  of  Leeway  and  wind 
angle  to  U  and  V  components  is  as  follows: 

U  Leeway  *  0.94  x  sin(163.6°)  «  0.266 

V  Leeway  *  0.94  x  cos( 163.6°)  *  -0.902 

The  leeway  and  current  are  applied  to  the  position 
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as  fol 1 ows : 


U 

V 

Leeway 

.266 

-.902 

Current 

.  » 756 

.JQZ 

Total 

1.032 

-.200  knots 

These 

values  are 

then  divided  by  the  distance 

between  gridpoints,  at  this 

latitude  ■  235.37  miles: 

.00439 

-.00085 

♦  position 

38.QQ787 

19.01413 

New  Position 

38.01226 

19.01328  for  1000Z  22  June. 

This  process  Is  repeated  hour  by  hour  until  the  last 
datum  has  been  calculated/  drift  plotting  ceases,  and 
calculations  for  position  error  and  Search  Radius  begin. 

6.  Position  Error 

Just  as  In  the  manual  method,  the  Initial 

position  error  of  the  survivor  craft  Is  computed  as  follows: 

x  ■  distance  of  Initial  (or  distress)  position  from 
last  known  position. 

A  ■  error  of  last  known  position 

B  ■  navigation  error  of  survivor  craft  as  a  percent 
of  distance  traveled 

The  Initial  position  error  ,  X,  Is 
X  •  A  ♦  B  x 

Since  Y  Is  the  nav'gatlon  error  of  the  searcher, 
and  de  is  the  drift  error,  the  Total  Probable  Error,  c,  of 
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Position  Is 


c  •  ^X*  ♦  V*  ♦  do* 

7.  SaftLX  flfilflf 

c,  Total  Probable  Error  of  Position  Is  xultlplled  by 
the  safety  factor  to  find  the  Search  Radius  as  prescribed  In 
the  Search  and  Rescue  Manual. 

I.  Sweep  Width 

A  numerical  approx  I  mat  Ion  to  the  date  In 

Richardson's  Sweep  width  tables  (Ref.  6)  nay  be  seen  to  be: 

W  •  -5.5*577  ♦  1.25117  log  (I)  ♦  2. 55**2  log  (¥)  ♦ 
.0102)  h 

where  I  Is  boat  length 

V  is  visibility  In  miles,  and 
h  Is  search  altitude  In  thousands  of  feet. 

In  addition,  corrections  to  Sweep  Width  for 
whltecaps  based  on  wind  speed  (v)  and  cloud  cover  (C)  can  be 
approximated  by: 

fl  •  .116  •  . 2961  -  .116  X*  (whltecaps,  X  •  v/10) 

F 2  •  1.1)1  -  . * 7 2 C  ♦  .1S2C1  (cloud  cover) 


These 

correction  factors 

are  multiplied 

by 

W 

to 

obtain  final 

Sweep 

Width.  This 

Sweep  Width  function 

Is 

applicable  to 

snal  1 

boats  only. 

No  Sweep  Width 

for 

large 

ships  or  for  llferafts  is  available  In  the  program. 


C.  VERSIONS  Of  THE  PROGRAM 

The  Search  Planning  Program  Is  written  In  FORTRAN  IV  for 
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us*  on  a  reaote  computer  terminal.  A  second  version  is 
written  for  use  on  the  operations  computer  at  Fleet 
Nuaerlcal  Weather  Facility.  Mo  reaote  terminal  operation  Is 
available  there,  so  these  two  versions  are  different  In 
several  ways: 

1.  I  mm t -Out nut 

a.  The  tla*  sharing  version  Is  written  for  entry  of 
data  by  the  search  planner  onto  a  reaote  console,  and  the 
results  are  printed  out  laaedlately  on  the  sea*  console. 
This  peralts  the  planner  to  control  the  program  during 
execution,  and  to  make  Intermediate  decisions  concerning  the 
problem.  An  eiample  of  the  output  of  the  time-sharing 
version  Is  shown  on  th*  neat  two  pages.  Note  that  I  teas 
typed  In  by  th*  user  are  shown  In  tower  case,  and  the 
replies  typed  on  the  console  by  the  computer  are  In  upper 
case.  The  final  output  Includes  sweep  width,  total  probable 
error  of  position,  search  radius,  and  calculated  wind. 

b.  The  batch* processing  version  at  FNWC  differs 
froa  this  In  that  It  reeulres  cards  to  be  punched  for  Input, 
and  thp  output  Is  printed  In  message  fora,  ready  for 
transailsslen  bach  to  tbe  reeuestor. 

For  the  operational  version  of  the  program  at  Fleet 
Numerical  Weather  Facility  the  following  Items  of  output 
shown  on  the  sample  tlm**sharlng  printout  were  not  mad* 
aval lablei 

lb 


# 


1  tar 

rirruTiov  orc.rs... 

MOMTTRT Y  s TARCm  RtANfll»‘C  PROrr**. 
fNTfB  irfwAY  copt,  1  RiAwr,  Turn  pr  sr  p  i  oj  |  o** 

>  1 1  1  I fnraf I 

mow  ma**y  START  i  wr  points  will  »*'  rvyrprri 

>  1 

l  I  ST  CMBOfOl  OC  |  r  A I  I  Y  :  ST*RTl»T  lints  •*•»'  ro<;  |  T  IP**S,  r»: 
"7<mn7  HAP,<A??.S*J,^P*OA.SW 

>  770900*  |un, 6067. On, 05715.  Pw 

60  67.0  S?  15. o  rwrrr;.'  * 

CMfOK.  TVPT  Y  PR  *' 

>  V 

fNTfB  HUMRfR  Pf  T I  Mf  S  PA  Tl»*  I S  TP  «T  CP*'PtrTTP  HOW: 

)  ? 

fNTfR  T*«fsr  DTPS .  r«:0?O9PO  MAR, 051050  MAP 

>  750100  Jun, 750700  Jun 

STARCH  1,  TRACK  l#  DATUM;  60  SO. 71  SI  S 1 . 96 

SfANCM  ?,  TRACK  J,  DATUM;  60  S5.76  SI  66.06 

fNTfR  LAST  KNOWN  POSITIOf.  r«:  55  7  7  .  ON,  0  76  S  5  . 0W 

>  6J00.  On, 05600.0** 

IS  l Iff  RAPT 

A  (l)BOAT,  (71SMIP,  OR  ( 5 ) A  I RCPAf T »( f NT f R  l,?,ORS> 

)  1 

f NTf R  SfARCw  UNIT  Tl*  f RROR ,  V.  f I :  S.O 
)  5.0 

I  NO  I  CATf  0  OISTPrSS  CRAfT  POSITION  TRPOr,!; 

I*  15.00(1)  •  0.15(7)  T  I  **r  S  P  I  $  T  .  f  RP**  IASI  KNOW*  "OS  I  T  I  ON . 
ANY  CORRfCTlOHS’*  (  T  Y  PC  0,1,7,  OR  1?) 

>  1 

rNTfR  PAPAMr Tf P( n : 

%  70.0 

IS  f  ACm  n  ATI**  Tl'*f  A  o»rrrprr*T  sr  ARP“  (  1 .  7 ,  5 , T TC  .  )  ? 

>  no 

fNTfR  SfARCM  ♦!!»♦*•»■*  roPorsPO»‘P  |  *>S  TO  PM  IP*  Tl**r: 

•  •  l  0  ST ARCw : 

>  1 

••70  STARCH: 


75 


i 


>  0 

TMF  rntiO»Ji»T.  A#f  TMF  ,'*Tir«  r>0 1 *' TS  rr>r  T«r 
1:  0  S FARC*: 

AB-50. 2H  51-57. OW  BOS.  FRBOB:  Al.t,  BAD  |  tts :  A*. 

W|Hf>  1%  .  ARTS  F«OM  Jf,  7  . 

SWFFB  WIDTH?  Y  OB  H 

>  Y 

FNTFR  SFARfM  ALT  I Tt  DF  (ri:IOO.) 

>  1000. 

fNTfR  VISIBILITY  IN  m|IFS  (FX:10.0) 

>  20.0 

CNTFB  BOAT  LENGTH  IN  FFFT  (fI:S6.0) 

>  10. 

FNTFR  CLOUD  COVER  IN  FRACTION  (FI:  0.0) 

>  0.1 

SWFFB  WIDTH.  A. 6  Mil rS. 

TMf  FOLLOWING  ART  TMf  DATUM  BOINTS  FOR  TMF 
7:  0  SFARCM: 

AO-55. 3N  51-AA.ftW  BOS.  FRROR:  Al.B,  RADIUS:  0 . 

WIND  3% . ARTS  from  26?. 

GOOD  |l»r«  ON  YOt't  SFARDMI  |  FT  Hr  cnow  |F  YOM  F|»*n 
THFLIFFRAFT 


IHC002I  STUB  fl 

R;  T-1.00/2.R7  13.ll.iR 


Hot* :  Input*  with  **>" 


HltFS 


Ml  LCS 


>1 


a.  Printout  of  Intermediate  track  positions  (SEARCH 
1,  TRACK  1,  etc.)  Is  not  produced. 

b.  On  scene  Wind 

Since  the  wind  printout  as  shown  In  the  sample 
t Ime-shar | ng  output  Is  only  the  wind  velocity  and  direction 
at  the  beginning  of  the  search,  and  may  not  be 
representat I ve  of  the  average  weather  conditions  throughout 
the  search,  this  value  Is  not  Included. 

c.  Sweep  width 

Since  the  the  computation  of  sweep  width 
requires  data  which  may  not  be  Invaedlately  available  to  the 
search  planner  at  the  time  of  his  request,  and  since  the 
formula  applies  only  to  small  boats,  this  feature  was  not 
Included  in  the  fleet  Numerical  version. 

d.  leeway 

for  the  fleet  Nim»erlcel  version,  leeway  Is 
calculated  as  described  in  D.2.  The  time  sharing  version 
uses  a  square  root  function  for  all  leeway.  If  N  Is  the 
percentage  drift  specified  (II  to  9t>,  the  leeway  function 
In  the  time  sharing  version  Is 

*tm4°  m 

This  function  produces  leeway  which  Is  equal  to 
N|  of  the  wind  for  wind  speeds  of  lt.O  knots.  for  other 
wind  velocities  the  function  follows  the  general  contour  of 
the  llferaft  leeway  (without  drogue)  table  In  the  Search  and 
Rescue  Manual . 
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III.  SEARCH  THEORY 


In  part  III  the  method  of  estimating  the  position  of  the 
search  object  and  the  rules  for  determining  the  site  of  the 
search  area  were  discussed.  This  Is  only  part  of  the  search 
planning  process;  this  section  will  consider  how  resources 
are  assigned  to  the  search,  and  how  the  effectiveness  of  the 
search  Is  measured. 

In  a  search  pettern  employing  parallel  sweeps,  the 
distance  between  these  sweeps  Is  called  track  spacing,  S. 
The  Search  and  Rescue  Manual  recommends  that  for  urgent 
caeet  the  track  spacing  used  be  eeuel  to  the  Sweep  width,  w, 
and  for  less  urgent  cases  It  be  equal  to  twice  the  Sweep 
Width.  Once  the  value  of  track  spacing  has  been  decided 
upon,  the  amount  of  searching  required  (In  miles  flown)  l 
can  be  computed  from  S  and  the  area  A  to  be  searched: 

l  •  A/S 

The  measure  of  effectiveness  of  a  search  which  Is  given 
In  the  Search  and  Rescue  Manual  Is  called  Probability  of 
Oetectlon,  and  It  Is  presented  as  a  series  of  curves  plotted 
agelntt  "coverage  factor"  or  W/S.  These  curves  are  called 
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“first  starch",  "second  starch",  etc.,  probabi 1 1 1 1 ts.  If  p 
is  tht  probability  of  dtttction  in  the  first  starch  curve, 
tht  probabi 1 1 1 Its  2,..n..S  art  given  by 

m 

p  •  1  -  (  1  -  p) 

Thtst  curves  do  not  take  Into  account  the  site  of  the 
area  starched,  however,  so  that  a  starch  can  be  conducted  In 
a  very  small  area  using  a  high  coverage  factor  which  will 
yield  a  very  high  probability  of  detection  according  to  the 
graph.  Common  sense  Indicates  that  to  confine  the  search  to 
a  very  small  area  will  probably  not  result  In  the  detection 
of  the  target.  Nowhere  does  the  Search  and  Rescue  Manual 
define  the  probability  that  the  target  Is  In  the  search 
area. 

The  absence  of  the  treatment  of  this  seemingly  critical 
factor  led  the  author  to  a  review  of  the  field  of  search 
theory  as  It  applies  to  survivor  search. 

The  basis  for  most  of  the  search  planning  criteria  In 
the  National  Search  and  Rescue  Manual  was  found  In  koopman't 
<19kk)  work.  Search  and  Screen! n*  (Ref.  2).  This  treatise 
developed  the  concept  of  sweep  width,  and  found  the 
probability  of  detection  for  various  situations. 

A.  Of TfCTION  LAMS 

.  Instantaneous  Probability  of  Detection 
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If  Y  is  defined  so  that  Che  probability  of  detecting 
a  target  In  a  short  period  of  tine  dt  is  Ydt,  then  the 
quantity  Y  is  called  the  I  nstantaneous  probability  density 
of  detection.  When  the  searching  is  done  during  a  tine  t 
under  constant  conditions,  the  probability  p(t)  of  detection 
is  given  by 

-  rc 

fit)  *  /  -  « 

(k) 

Search  and  Screening  developed  this  quantity  for 
various  situations. 

a.  Definite  Range  Law 

This  law  assumes  that  there  Is  a  certain  range  R 
from  the  observer  beyond  which  It  Is  Impossible  to  detect 
the  target;  Inside  this  range  detection  Is  certain. 
Instantaneous  probability  of  detection  by  this  law  Is 

>  tr)  -  ou  rc  ^ 

_  (S) 

X  lr)  .  0  r  >  R 


b.  Inverse  Cube  Lew 

The  Inverse  cube  law  Is  the  result  of  making  the 
following  assiaeptlons: 

(1)  The  observer  Is  at  an  altitude  h  above  the 


oceen. 

Its  wake. 


(2)  The  observer  dt ‘acts  the  target  by  seeing 


(S)  The  Instantaneous  probability  of  detection 


Is  proportional  to  the  solid  angle  subtended  at  the  point  of 
observation  by  the  wake. 


The  calculation  of  the  solid  angle  is  shown  In 

figure  2. 


Figure  2.  Calculation  of  Solid  Angle 


Solid  Ang I e  ■  >>•  /3 

ot-  C/s  ,  0  ’  ^  h 

but  cA’v  *  ^/5 

thus  C  *  *  /*  - 

Solid  Angle  -  tilth/*3  -  a>h/(r2«-h  ) 

Since  y  was  assumed  to  be  proportional  to  the 

sol  id  angle,  ,  . 

*  r 

where  k  included  all  factors  such  as 
meteorological  visibility,  etc. 
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Assuming  further  that  in  the  majority  of  cases  r 


is  much  larger  than  h,  then 

<  * 

V  -  — 

2.  Lateral  Range 

When  the  observer  and  the  target  are  on  their 
straight  courses  at  constant  speeds  for  a  long  time  before 
and  after  their  closest  point  of  approach  the  probability  of 
detection  may  be  shown  to  be  a  function  only  of  the  distance 
at  closest  point  of  approach.  This  distance  is  called  the 
lateral  ranee,  and  is  denoted  as  x.  This  probability  of 
detection,  p(x),  plotted  against  lateral  range,  is  called 
the  lateral  ranee  curve. 

Twice  the  area  under  the  lateral  range  curve  Is 
defined  as  effective  search  width  (or  Sweep  Width), 


r  Cn 


Koopman  shows  that: 


(6) 


"The  effective  search  width  is  twice  the 
range  of  a  definite  range  law  of  detection 
which  Is  equivalent  to  the  given  law  of 
detection  in  the  sense  that  each  of  the  two 
laws  detects  the  same  number  of  uniformly 
distributed  targets  of  Identical  velocity." 

Thus,  for  the  definite  range  law  of  detection  p(x)  ■ 
1  If  xi  R,  and  p(x)  »  0  If  x  >  R,  so  that  W  ■  2R.  To  find  W 
for  the  Inverse  cube  law,  first  p(x)  may  be  shown  to  be 


p(x)  =  |  -  c 


(7) 


2  k  h 
co1- 


where  w  Is  searcher's  speed. 


To  find  the  probability  of  detecting  a  target  known 
to  be  In  an  area  A,  where  both  searcher  and  target  are 
moving  over  the  ocean  In  compl I cated  paths,  the  situation 
called  random  4£a££li  may  be  considered.  This  requires  the 
following  assumptions: 


a.  The  target's  position  Is  uniformly  distributed 
In  A,  given  it  has  not  been  detected 

b.  The  observer's  path  Is  random  In  A  In  the  sense 
that  It  can  be  thought  of  as  having  Its  different  (not  too 
near)  portions  placed  Independently  of  one  another  In  A. 

c.  W  Is  small  with  respect  to  the  area  A. 

If  the  searcher's  path,  of  length  l,  can  be 
considered  to  be  composed  of  n  separate  and  equal  portions 
of  length  L/n,  then  the  probability  that  the  target  be 
detected  in  any  one  path  Is  Wl/nA.  The  probability  that  the 
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target  not  be  detected  at  all,  then  is 


m.  \r 
/ 


and  for  large  n, 

-vvl/A 

p 1  I  -  t 

This  Is  the  formula  of  random  search. 


( 10) 


(11) 


The  most  common  method  of  non-random  searching  is 
one  which  uses  a  series  of  straight  parallel  sweeps  a 
distance  S  (or  track  spacing)  apart,  where  the  position  of 
the  target  is  fixed  with  relation  to  the  sweeps.  Let  the 
target  be  uniformly  distributed  in  an  area  of  size  A. 
Search  and  Screening  gives  the  following  result  for  P(L), 
the  probability  of  detecting  with  path  length  L*A/S. 

(1)  Definite  Range  Law 

Since  detection  will  surely  occur  if  the  target 
Is  within  the  detection  range  W/2,  then  the  probability  of 
detection  is  equal  to  the  coverage  factor  W/S. 


Thus, 


=  wl/a 


L_  1  A/w 


L  >  A/W 


(12) 


(2)  Inverse  Cube  Law 


„  (-f—)  -  erf  U  A  I 


(13) 


(3)  Uniform  Kindon  Search 


From  (  11 ), 

f  =  ‘  -  * 

or  when  S*L/A, 

-  wv  /  i 

?  •  -  *  (1%) 
Tht  thro*  probabilities  using  parallel  sweeps 
are  shown  In  Figure  3. 


Figure  3.  Probability  of  Detection 

I.  D I  ST*  I  BUT  ION  OF  SEARCH  EFFORT 

1.  fitfllfAl 

Another  result  In  Search  and  Screening  Involves  the 
definition  of  search  effort  density  as  a  function  of  («,y) 
In  A. 


If  •  Is  the  area  of  a  subregion  in  A,  and  l  Is  the 


length  c'  observer's  path  in  B,  then  the  eap'ession  L  /B 


is  defined  to  be  the  observer's  track  density  in  B.  As  B 
becomes  very  smell  (approaches  a  point  (a#y))  the  search 
dens  1 1  y  {  (  * ,  y )  is  the  limit  of  ML  /B,  or  Ml(a,y)#  where 


(  IS) 


//  V  <v  y> 


d  a  <i  y  -  VV  t.  -  £ 


(  It) 


Thus  If  (a,y)  is  applied  at  (i,y),  the  probability 


of  detecting  a  target  located  there  is 

-  <e(*ty) 

I  _  i 

P 


I  -  C. 


(17) 


Given  the  search  density  function  ^'(a.y)  and  the 
terget's  distribution  p(a,y)#  the  prooabillty  of  detecting 


the  target  Is  -  4  (0 1  y) 

r(r)  -  ((rt'.th  1  —  e 


(  II) 


The  ^  (a,y)  which  maalmlaes  f(  </  )  is  shown  by 
koopman  to  be 

r(’,  i) 


(f  (*  ,y )  J  tog 


( if) 


when  (a#y)  Is  In  A,  and 

f  (*,  y )  •  °  (20) 

when  (  a  #  y )  Is  not  In  A  , 

where  A  denotes  that  portion  of  A  In  which  the 
eapresslon  (  >9  )  Is  non-negat I ve.  The  value  of  b  is 
determined  by  using  equation  (II). 


7.  Hormal  Tifttt  Distribution 


If  the  target's  distribution  Is  a  bivariate  circular 
normal  distribution  about  a  point (0,0)  In  a  cartesian 
coordinate  system, 

P1’-)1)  -  z7o">  t 

(Jl) 


r  ‘ 


Is 


then  the  search  effort  density  which  maximizes  P( f  ) 

‘  r  <  u. 

(22) 

and 

(f('.y)  *  0  r  >«- 

on 


Thus  letting  a  be  the  radius 
assumes  the  critical  value  b,  k 

A  ir  *  tt  iw  ,  b  5  nrs  • c  * 

where 


*  - 


[•*<’•  £l  4 

tt 


r  at  which  p(  i,  y ) 

(2%) 

(25) 


The  form  of  this  equation  Is  an  Inverted  parabola  of 
revolution  with  Its  center  at  (0,0)  with  a  maximum  value  of 
*‘/2dl. 


57 


IV.  TMf  UNI*0P*4  COVERAGE  SQUARE  SEARCH 

Although  (  19  )  gives  the  ’’optimum**  search  effort 
distribution,  this  optimum  distribution  Is  difficult  to 
achieve  In  practice.  flying  In  a  large  circular  area  with  a 
continuously  varying  effort  is  no  simple  navigational  task. 
Much  easier  Is  the  common  practice  of  (1)  defining  a  square 
area  with  Its  center  at  the  center  of  the  target's 
distribution  and  (?)  conducting  a  regular  search  of  uniform 
density  throughout  this  area.  If  the  target  is  not  found, 
then  (3)  the  area  Is  expanded  and  the  search  repeated. 
Although  a  single  search  of  this  type  does  not  conform  to 
the  parabolic  optimum  search  effort  Ulstr Ibut Ion,  a  series 
of  them,  each  larger  than  the  one  before,  may  be  a  very  good 
approalmatlon. 

The  most  Important  advantage  of  this  plan  Is  Its  ease 
of  navigation,  allowing  the  pllot(s)  free  time  to  look  for 
the  target. 

A.  THE  PROBLEM 

The  problem  was  to  determine  the  optimum  site  of  a 
square  area  within  which  a  uniform  search  with  parallel 


sweeps  Is  to  be  conducted 


ft.  TERNS 


Th«  i «rm  and  definitions  in  the  Search  and  Rescue 
Manual  differ  from  those  In  Search  and  Screening.  In  the 
Search  and  Rescue  Manual,  probability  of  detection  is  used 
only  for  the  probability  of  detection  with  parallel  sweeps, 
(given  a  uniform  target  distribution  within  the  search 
area).  In  order  to  avoid  confusion,  the  following  terms  are 
defined: 

P  -  Probability  of  Success  -  The  probability  that  on  a 

given  search  the  target  will 
be  found. 

p  -  Location  Probability  -  The  probability  that  a  target 

Is  within  a  specified  area. 

d  -  Probability  of  Detection  -  The  probability  of 

detecting  a  uniformly 

distributed  target  given 
parallel  sweeps. 

w  -  Sweep  Width  -  The  effective  search  width  of  the 

searcher;  the  area  under  the 
lateral  range  curve. 

S  -  Track  Spacing  -  The  mean  distance  between  adjacent 

parallel  sweeps  of  the 

searcher. 

I  -  Path  length  -  The  total  number  of  miles  flown  by  the 

searcher  within  the  search 

area. 

c  -  Probable  Error  of  Position  -  Probable  error  of  the 

target's  position  from  datum; 
given  a  normal  dlstr (button, 
c«l.ll  ■  standard  deviation  of 
the  target's  distribution. 

R  *  Search  Radius  -  The  search  radius  Is  defined  as 

one-half  the  side  of  a  square 
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search  area. 


C.  THE  SOLUTION 


1.  IlMLJUaifcilH 


The  objective  is  to  maalmlxe  the  Probability  of 
Success  P,  where  the  detection  function  for  a  uniformly 
distributed  target  known  to  be  in  the  area  Is  defined  as 
d(M,l,R,a,y).  P  Is  expressed  as 


For  a  uniform  search,  the  probability  of  detection  is 
constant  and  equal  to  d(w,l,R)  Inside  the  search  area  A,  and 
Is  taro  outside,  (given  uniform  visibility  and  other 
factors),  and  so  d(W,l,R)  can  come  out  of  the  Integral: 


It  Is  convenient  to  choose  a  and  y  as  distances  In  a 
cartesian  coordinate  system  with  the  target's  most  likely 
position  (datum)  at  (0,0). 


2.  Calculation  of  Individual  Tarm* 

Based  on  the  Inverse  cube  law  for  parallel  sweeps. 


kO 


i. 


toopiMn  dtrlvts  tht  following  piprpttlon  for  probability  of 


dotPCtion,  fro«  ( • 3  > : 


* 


) 


whprp  S  •  A/ 1  •  (  ?R  J*/l. 
TM»  It  ooulvalpnt  to  thp  M 
dp t PC  1 1  on  glvpn  In  thp  Nptlonpl 
Tbp  location  probability 
for  tbp  circular  blvarlata  norma 
//  » 


flrtt  Saarcb"  probability  of 
Saarcb  pod  Ratcua  Manual. 

,  P.  It  glvpn  by  tbp  formula 
I  dlttr  Ibutlon: 

iv  * 

cl  i  y 


Whprp  R  dpflnpt  P  tQuPrp  cpntprpd  pt  (0,0)  with  tldpt 


flgurp  b  it  a  graph  of  tbit  function,  along  with  * 
tlmllar  cur vp  t bowing  tbp  cumulatlva  probability  for  a 
clrcla  wltb  radlut  R.  Slnca  tbp  touarp  arpa  It  actually 
largpr  (It  It  c  I  renter  I  bad  about  tbp  clrcla)  Itt 
probability  for  any  radlut  R  It  groator.  Howavpr,  If  tbpta 
cunwjlatlvp  probabl  I  1 1  lot  arp  plot  tad  agalntt  tbp  araa  of 
tbp  taarcb,  tbp  probability  for  a  touarp  will  ba  tllgbtly 
latt  than  for  a  clrcla  (tap  flgurp  S).  Although  thorp  It 
tomp  lott  of  pfflclpncy  In  taarchlng  a  touarp  Inataad  of  a 
clrcla,  thlt  lott  It  vary  Minimal,  about  It. 


S«w«r«  Plotted  ago lost  Area. 


NOTE:  Although  the  standard  deviation  is  the  commonly 

accepted  measure  of  variability,  the  National  Search  and 
Rescue  Manual  uses  the  term,  "Total  Probable  Error  of 
Position",  Its  value  Is  computed  using  estimates  for  drift 
error,  target's  position  error,  and  search  craft  navigation 
error.  These  errors  are  squared  and  added,  and  then  a 
square  root  Is  taken  to  produce  Total  Probable  Error.  This 
term  Is  not  specifically  defined  in  that  publication,  so 
that  one  of  the  objectives  of  this  solution  was  to 
demonstrate  that  It  is  equivalent  to  the  statistically 
well-defined  term  'Probable  Error',  which  is  a  median  error: 
that  is,  half  of  all  errors  have  absolute  values  greater 
than  that  error,  and  half  have  absolute  values  which  are 
smaller.  The  absence  of  this  definition  is  a  gap  in  the 
present  search  planning  process,  because  planners  are  called 
upon  to  make  estimates  of  errors  without  knowing  how  these 
errors  are  defined. 

>'  I  Hi  fiiOTil  Sal  ml  an 

To  determine  the  optimum  site  for  a  square  search, 
given  a  certain  amount  of  searching  miles  available  (l), 
sweep  width  (M),  and  probability  distribution  (p(a,y)),  a 
straightforward  maalmliatlon  method  Is  employed  In  the 
following  sentience: 

a.  Determine  total  effort,  E. 


b.  Choose  some  arbitrary  initial  search  radius,  x. 

c.  Determine  the  search  area.  A,  for  this  x: 

A  ■  4  x* 

d.  Determine  effort  density  for  a  uniform  search: 


e.  Determine  probability  of  detection  d  from  (3) , 

f.  Determine  the  probability  p  from  (2>;  with  R  * 


g.  Determine  the  probability  of  success, 

P(x)  ■  p( x)d( x) . 

h.  Increase  x  by  Ax  and  repeat  calculations.  If 
P(x)  increases,  continue  adding  Ax  until  a  local  maximum  Is 
encountered.  If  P(x)  decreases  Initially,  reduce  x  by  2Ax 


and  continue 


reached. 


Ax  until  a  local  maximum  is 


Since  a  cumulative  probability  distribution  Is 
continuous  and  has  a  slope  throughout  its  range  of  2.  and 
since  the  probability  of  detection  function  d(x)  Is 
continuous  with  a  slope  of  <  0,  the  function  P(x)  »  p(x)d(x) 
can  be  shown  to  be  unlmodal,  and  so  a  local  maximum  in  P(x) 
will  be  an  absolute  maximum. 

4.  Add! lcatlon 

The  first  use  of  this  method  was  to  check  the 
National  Search  and  Rescue  Manual's  recommended  search 
planning  guidelines  to  see  if  they  produce  square  searches 
of  optimum  sizes.  These  guidelines  are  given  in  the  form  of 


) 


'Safety  Factors',  or  search  radii,  for  repeated  searches,  as 
multiples  of  Total  Probable  Error  of  Position,  e.g.. 


Search 

1 

2 

3 

4 

5 


Safety  Factor 
1.1 
1.6 
2.0 
2.3 
2.5 


As  far  as  effort  density  is  concerned,  there  are  two 
recommendations  given  in  the  SAR  Manual.  One  is  for  a  track 
spacing  equal  to  sweep  width,  for  urgent  cases,  and  the 
other  for  a  track  spacing  equal  to  twice  the  sweep  width, 
for  less  urgent  searches.  Both  of  these  recommendations 
were  tested. 

The  approach  used  was: 

a.  Set  some  arbitrary  value  for  W 

b.  Set  S*W  so  that  WL/A.  -1. 

c.  Find  R  from  "safety  factor"  table  for  the 
appropriate  search  (starting  at  1). 

d.  Find  path  length,  L: 


e.  Beginning  at  x  ■  R,  and  with  W  and  L  above,  use 
the  general  solution  process  described  In  paragraph  3.  to 
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the  optimum  search  radius  for  that  path  length 


f i nd  x ' , 

f.  Assuming  that  the  optimal  search  Is  executed 
unsuccessfully,  reevaluate  the  target's  distribution  p(x) 
(see  Para.  5.)  and  repeat  the  process  from  c.  above.  After 
five  successive  searches  have  been  computed,  stop  the 
computation. 

g.  Compare  the  resulting  five  optimum  sizes  with 
the  sizes  recommended  in  the  Safety  Factor  table  from  the 
Search  and  Rescue  Manual. 

5.  Bavesian  Analysis 

Bayes'  Law  of  Probability  says  that  an  estimate  of 
the  probability  of  a  hypothesis  can  and  should  be  revised  on 
the  basis  of  the  outcome  of  a  related  event.  In  the  case  of 
search  theory,  the  hypothesis  Is  that  the  target  Is  Inside  a 
given  area.  A  search  is  conducted,  and  the  event  Is  either 
the  detection  or  the  non-detection  of  the  target.  This 
event  should  revise  the  initial  estimate  of  the  probability 
that  the  target  was  In  the  area.  Of  course.  If  the  event  is 
detection,  the  new  probability  that  the  target  is  in  the 
area  Is  1.  if  the  event  is  non-detection,  then  Bayes* 
formula  must  be  used  to  calculate  the  new  probability. 

1.  Example 

The  probability  that  a  target  is  in  a  certain  search 
area  has  been  determined  to  be  50$,  and  a  search  is 
conducted  In  this  area  with  a  uniform  probability  of 
detection  of  50$.  The  resulting  probability  cf  success  is  .5 
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x  .5  «  .25  or  25:. 

For  the  use  of  Bayes*  formula,  the  following  terms 


are  defined: 


A 

The 

hypothesis  that  the 

area . 

target  is 

In  the 

A 

The 

hypothesis  that  the 
the  area. 

target  is 

n^t  in 

0 

The 

event  of  detectior. 

D 

The 

event  of  non-detect I  on 

• 

Then 


p(A  ID)  = 


v  IpU) 

P^lA)  v  p  ( A  )  •  p(D|A) 


When  the  target  Is  detected: 


C. 5"  *  C  h  |  o 

OS  *  0-5T  r  C.‘j  *  o.  r> 


as  expected.  If  the  target  Is  not  detected,  the  new 
probability  would  be 

p(a) •  p(d\ a) _ 

p(A) '  r(5  M)  f  f(a  )•  p(  5 1  a  ) 


os  »  os _  ! 

C*  S'  ,  0-5  T  O.X  A  I.  o  ~  ~3 
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Whereas  the  target  had  a  50$  a  priori  probability  of 


being  In  the  search  area.  Its  new,  or  a  posteriori 
probability  of  being  in  the  area  is  33.3$. 

For  the  solution  in  paragraph  4.,  the  entire 
probability  distribution  of  the  search  object  was 
reevaluated  between  successive  searches.  This  was  done  for 
p(x)  inside  the  search  area  by 

»'(*)  -  - 

r  h  -  p(*)cO  <*«> 


and  for  p(x)  outside  the  search  area  by 

r'U)  -  l-f-U)  rp'W 


(27) 


where  R  is  the  radius  of  the  search  conducted. 


in  order  to  illustrate  the  general  form  of  the 
probability  distributions  between  searches,  a  function  was 
used  which  is  similar  to  a  density  function.  This  function, 
f(x)  was  defined  to  be 


f  Cx) 


/  d(p(t)) 
\  cl* 


(28) 


where  x  is  the  measure  of  radius,  or  again,  one-half 
of  one  side  of  a  square,  and  8xdx  is  the  differential  of 
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area  (Figure  6  ) . 


8  “  ^  A. 


Figure  6.  Derivation  of  Density  Function  for  Square 

Area 

Curve  1  in  Figure  7  is  the  density  function  before 
any  searching  is  done,  as  derived  by  (28)  from  (21).  Curves 
2  through  5  are  these  functions  after  successive  searches. 

Figure  8  is  a  similar  set  of  curves  showing  the 
cumulative  probability  distributions. 

6.  Resul ts 

This  method  (from  4.)  gives  the  following  optimum 

searches: 


Search 

1 

2 

3 


Optimum  Radius 
1.08 
1.64 
2.07 


Safety  Factor 
1.1 
1.6 
2.0 
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Probability  Density 


r 


e. 


\ 


\ 

\ 

\ 

\ 

\ 


\ 


R  / 


Figure  7.  Probability  Density  with  Application  of  Bayes'  Law 
Between  Successive  Optimum  Square  Searches. 
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Figure  8.  Cumulative  Probability  with  Application  of  Bayes'  Law 
Between  Successive  Optimum  Square  Searches. 


4 


2.43 


2.3 


5 


2.74  2.5 


Since  these  results  were  very  close  to  the  'Safety 
Factors'  listed  above.  It  was  concluded  that  (1)  the  safety 
factors  were  oriflnally  food  calculations,  and  (2)  they  were 
calculated  usinf  the  mathematical  definition  of  probable 
error. 


It 

should  be 

noted,  however. 

that 

these 

safety 

factors 

are  close  to 

optimum  search 

radi  1 

only  for  the 

condl t Ion 

where  track 

spacing  equals 

sweep 

width. 

The 

second 

recommendation. 

track  spacing 

equa  1  s 

twice 

sweep 

width. 

does 

not  produce 

anywhere  near 

optimum 

search 

areas. 

When  a  track  spacing  S'  equal  to  twice  the  sweep  width  is 
used,  the  new  path  length  L'  is  given  by 

-  ±  L 

L  ”  5'  ;w  2.  u 
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When  there  Is  less  search  effort  available  the  size  of 
the  search  area  as  well  as  the  effort  density  should  be 
reduced.  The  SAR  Manual  fives  no  guidance  as  to  how  much  the 

area  should  be  reduced  with  a  given  amount  of  available 

effort.  Conversely,  there  Is  no  Indication  as  to  how  excess 
search  resources  should  be  used,  if  such  a  condition  should 
exist. 

What  Is  needed  Is  some  way  for  the  search  planner  to 
determine  the  optimum  search  configuration  directly  from  the 
amount  of  search  effort  available. 

D.  THE  GENERAL  CASE 

Assume  that  a  series  of  optimum  square  searches,  each 
larger  than  the  preceding  one.  Is  a  good  approximation  of  an 
optimum  search  effort  distribution.  Consider  all  the 

searches  of  the  series  combined  Into  one  search  effort,  so 

that  the  search  effort  density  0(r)  assumes  a  step  function 
which  is  close  to  the  Ideal  effort  density. 


The  radius  r  of  the  combined  search  Is  equal  to  the 


radius  of  the  last  (and  the  largest)  individual  search  in 
the  sequence. 

If  the  search  had  been  a  circular  one,  and  if  the  effort 
density  had  conformed  to  the  optimum,  the  expression  for 
this  radius  would  be 


a 


A  cr1  £ 

rr 


i 

4 


But  since  the  search  is  composed  of  a  small  number  of 
square  searches,  this  same  formula  is  not  expected  to  apply. 
There  should,  however,  be  some  way  to  approximate  a  function 
which  will  define  the  optimum  radius  based  on  available 
effort  and  the  target's  distribution.  The  simplest  way  to 
do  this  is  to  compute  total  search  effort  for  a  few  optimum 
square  searches,  plot  the  radius  of  each  of  the  searches 
against  the  cumulative  search  effort,  and  attempt  to  smooth 
a  curve  through  them. 

This  was  done  for  the  ten  searches  (five  at  W  *  S  and 
five  at  W  *  2S)  calculated  above,  and  the  resulting  curve  is 
shown  In  Figure  9,  Just  as  radius  Is  in  terms  of  probable 
error,  so  are  the  factors  W  and  L  of  total  effort 
'normalized'  by  dividing  each  by  probable  error.  Thus 
'normalized'  effort  is  a  multiple  of  the  square  of  probable 

It  should  be  cautioned  that  the  curve  is  only  an 
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approximatl on,  and  will  not  be  valid  if  the  effort  available 
for  a  single  search  is  very  large.  Thus  if  a  search  is 
computed  which  shows  a  coverage  factor  (W/S)  greater  than 
1.8,  then  the  results  of  the  curve  will  not  give  an  optimum 
search,  and  some  other  method  will  be  required.  Again,  the 
effort  plotted  is  cumulative,  so  that  all  the  effort  exerted 
in  earlier  searches  must  be  added  to  the  effort  available 
for  the  present  search  before  entering  the  table. 

Example:  One  150-knot  aircraft  is  available  for  six 
hours  of  searching;  the  total  probable  error  of  position  has 
been  calculated  to  be  31.9  miles,  and  W  is  estimated  at  10 
miles. 

Total  effort  then  is  found  as  follows: 

E»  150kts  X  6hrs  X  lOmiles(W) 

■  9000  square  mi les 
In  terms  of  probable  error,  (c»31.9) 

E( norma  1  I  zed )  ■  9000/(31.9  X  31.9)  -  8.9 


From  the  graph,  r  ■  1.4, 


R  ■  1.4  X  c  ■  44.6  miles,  or  89.2  miles  on  a  side,  and 


S-4XRXR/L  -  8.75  mi les. 
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If  the  target  is  not  found  on  the  first  search,  and  a 
second  is  to  be  conducted,  the  amount  of  available  effort  is 
calculated  for  the  second  search.  The  normalized  effort  for 
this  search  must  be  combined  with  all  previously  expended 
effort  before  entering  the  graph.  Thus  if  t ( norma  1 i zed )  for 
the  second  search  is  calculated  to  be  10.0,  the  graph  should 
be  entered  with  a  figure  of  10.0+8.9  or  18.9  to  find  the 
radius  of  the  second  search.  Each  search  area  will  be 
larger  than  the  preceding  one. 

This  decision  curve  should  prove  to  be  a  more  flexible 
tool  for  search  planning  than  the  Safety  Factor  table,  since 
it  provides  a  continuous  range  of  search  radii  for  different 
amounts  of  effort  available. 
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V.  THE  SECTOR  SEARCH 


There  Is  one 

search  pattern 

In  general 

use 

which  does 

effect  a  greater 

concent  rat  Ion  of 

effort  at 

the 

center  of 

the  search  area,  and  that  It  the  sector  search.  This  starch 
patttrn  Is  executed  by  a  series  of  sweeps  crossing  over  the 
center  of  the  search  area  at  various  angles.  These  sweeps 
are  connected  at  their  endpoints  by  chords  flown  near  the 
perimeter  of  the  area.  This  pattern  covers  a  circular  area 
and,  perfectly  navigated,  produces  a  track  as  shown  In 
Figure  10. 

A.  SEARCH  EFFORT  DISTRIBUTION 

The  total  effort  expended  In  a  sector  search  can  be 
expressed  by  the  Integral 

£  c  ^  <*>(')  2"  r  dr 

where  R  is  search  radius 

r  Is  the  distance  from  the  center  of  the 

search 
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(fir)  It  the  search  effort  density  at  some 

distance  r. 

This  density  function  (fir)  Is  simply  W  times  the 


number  of  peths  which 

Intersect 

a  circle 

of 

radius  r. 

mu  1 1 1  pi  led  by 

dr  to  give 

total  path  length 

between  r  and 

r«dr,  and  then 

divided  by 

the  area 

between  r 

and 

r*dr.  Thus 

(fir )  ’ 

Wn  dr 

}  IT  r  i‘r 

- 5 

Wn 
in  r 

(19) 

Then, 

E  1  WL  * 

I  W" 

/  itrr 

i  n  r  dr 

c* 

c 

11 

(30) 

'0 

and, 

L«nR 

as  expected. 

Here  n  Is  the  number  of  ' spokes'  produced,  or  twice  the 
number  of  complete  sweeps  across  the  circle.  This  density 
function  Is  Illustrated  in  Figure  11. 

B.  PATTERN  ERROR 

There  will  always  be  some  error  In  flying  any  search 
pattern,  which  should  be  taken  Into  account  for  the  sector 
search.  The  navigation  error  of  the  searcher  has  already 
been  taken  Into  account  when  computing  the  Total  Probable 
Error  of  Position.  The  pattern  error,  however,  can  be 
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Distance  from  Center,  R 

Figure  11.  Search  Effort  Density  in  e  Perfect 
Navigated  Sector  Search, 
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different  from  this  navigation  error .  Whereas  navigation 
error  contldart  only  the  searcher's  ability  to  find  the 
datum  point  In  order  to  daflna  tha  taarch  area#  pattarn 
error#  at  It  It  utad  here#  datcrlbat  hit  flying  error  within 
the  area.  If#  once  tha  taarchar  flndt  what  he  thlnkt  to  be 
datum#  he  drops  a  snoke  float  or  tome  other  marking  device# 
then  hit  pattern  error  will  be  much  imatler  than  hit 
navigation  error.  Conversely#  If  when  arriving  In  the 
search  area  ha  stops  referlng  to  his  navigation  aids  and 
relies  thenceforth  on  dead  reckoning  only#  his  pattarn  error 
will  be  greater  than  his  navigation  error. 

Although  by  Its  nature#  a  flying  error  It  bivariate 
(Ilka  the  error  of  target  position)  a  simplification  was 
made.  Since  the  pilot  Is  flying  long  sweeps  which  have 
their  endpoints  relatively  far  away  from  the  area  of  concern 
(the  center  of  the  area)#  only  his  lateral  position  error 
(which  Is  univariate)  was  considered. 

It  was  assumed  for  this  analysis  that  the  lateral 
component  of  pattern  error  Is  nearly  normal#  and  that  it  can 
somehow  be  calculated  or  estimated. 

The  expression  for  the  total  effort  expended  in  flying 
the  spokes  of  a  sector  search  Is 

e=  lWn//  /TT®,  £  rdrJ*  <51) 

//  (r>*n<x  .  leterel 

(  O' p  -■  O'  Lateral  p.Mtvro  trr(r  ) 
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The  derivative  of  this 
used  as  the  search  effort 
search  with  pattern  error. 


U  ('  i 


r,/> 

i 

Jzrr  ©7 


t 


expression  with  respect  to  r  was 
density  function  for  a  sector 


f 

t  i '  h  ‘  x 

ilc*  (32) 


For  the  analysis  this  basic  density  function  was 
computed  with  sweep  width,  standard  deviation  of  pattern 
error,  and  number  of  spokes  all  equal  to  1.  To  find  the 
density  function  for  any  particular  set  of  parameters,  the 
density  is  multiplied  by  W  and  n,  and  the  horizontal  scale 
is  multiplied  by  the  standard  deviation  of  the  pattern 
error. 


C.  OPTIMUM  SECTOR  SEARCH 

The  plot  of  this  density  function  is  found  in  Figure 
12.  An  inspection  of  this  plot  will  show  that  it  is  close  to 
a  quadratic  function  out  to  r  ■  3  standard  deviations. 
Since  Koopman's  Ideal  search  is  a  quadratic,  some  guidance 
was  provided  as  to  the  optimum  radius  for  the  sector  search, 
based  on  the  pattern  error.  It  was  then  necessary  to  adjust 
the  number  of  spokes  to  arrive  at  the  Ideal  effort 
expenditure  to  match  this  radius.  Koopman's  Ideal  search 
effort  density  Is  a  function  of  search  radius  and  target's 
position  error: 


6b 


Search  Effort  Density 


R/Pattern  Error  in  o- 

Figure  12.  Search  Effort  Density  in  a  Sector  Search 
with  Pattern  Error. 
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( 

and  the  radius  Is  a 


and  available  effort: 

jc ' 1  EL 


a  ; 


if 


1. 


t  *£•  1 1 

t  >  <*■ 

function  of  target  position  error 


Solving  for  effort,  E: 

TT 

E  = 


W  L 


However,  In  the  case  of  the  sector  search,  L,  the 
number  of  miles  flown.  Is  simply  the  number  of  spr kes,  n, 
times  the  radius.  The  chords  flown  at  the  perimeter  are 
disregarded  In  this  analysis: 


L  ■  a  X  n 


Substituting,  and  solving  for  n: 

0?  TT 

h  -  - 

4crz  W 

Ideally,  then,  since  the  optimum  a  Is  known  to  be 
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Table  1  is  a  tabulation  of  this  function  for  various 
conditions  of  total  probable  error  of  positloi.  c,  and 
pattern  error.  The  n  derived  from  the  table  must  be  divided 
by  sweep  width  W  lg  obtain  number  of  ipokes. 

From  the  example  of  the  sunken  fishing  boat#  c  ■  31,9 
and  W  ■  10.0  miles,.  If  the  searcher's  pattern  error  were  .’3 
miles,  then  search  radius  P<  would  equal  58,  5  miles. 

Interpolating  in  the  table,  n  ■*  100,  then 
Number  of  spokes  **  100/10  •  10 

and  the  angle  between  adjacent  search  legs  *  3E0/10  *  36 
degrees . 

Total  search  miles  flown  would  be 
10  X  58,5  -  585 

plus  the  amount  of  flying  at  the  perimeter  (in  order  to 

fly  me  pattern,  the  aircraft  must  traverse  a  distance 

somewhat  less  than  one-half  the  circumference  of  the  area,': 
585  ♦  58,5  X  IT  >  53  5  ♦  li’4  •  7  70  miles. 

This  table  for  sector  search  will  not  be  useful  until  an 
accurate  means  of  pr*d  i  r.  t.  \  *,g  pattern  error  is  found.  The 

function  ^iven  >  « 3d  )  Is  extremely  sensitive  to  pattern 


erroi,  thus 

accuracy 

for  this 

parameter 

Is  required 

to 

use  the  table. 

Even  If  It 

1  s  not 

actua 1 1 y 

used  to 

determl ne 

the 

parameters  of  a 

sector 

search. 

the  table 

does  serve 

to 
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Table  1.  Optimum  Sector  Search.  Oivlde  quantity  In  "N  x  W"  table  by  Sweep 
Width  to  obtain  number  of  search  lees. 


i I  I u * t r a t •  that  accurate  navigation  does  not  necessarily 
product  the  optimum  sector  search,  and  that  onty  under  the 
unusual  conditions  of  an  extremely  accurate  estimate  of 
target  position  Is  It  deslreable  to  reduce  pattern  error. 
On  the  contrary.  In  most  Instances,  effectiveness  may  be 
enhanced  by  Incraai I  ns  pattern  error. 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  THE  SEARCH  PLANNING  PROGRAM 

It  Is  recommended  that  the  Coast  Guard  begin  Immediately 
to  make  full  use  of  the  Search  Planning  Program  now 
operational  and  undergoing  testing  at  the  Fleet  Numerical 
Weather  Central ,  Monterey.  A  guide  for  users  of  this  program 
is  included  as  Appendix  A. 

1.  Precautions 

The  following  cautions  are  appropriate  for  users  of 
the  Information  provided  by  this  program. 

a.  Region  of  Validity. 

Because  the  weather  data  points  used  In  this 
program  are  approximately  200  miles  apart,  drift 
calculations  provided  for  areas  within  100  miles  of  shore 
should  be  doubl e-checked  by  hand  before  using.  If  hand 
calculations  produce  a  larger  value  for  total  drift,  then 
the  hand-calculated  answers  should  be  used  in  lieu  of  the 
computer  solution. 

b.  Time  of  Datum 

Weather  predictions  in  the  format  required  for 
this  program  are  not  produced  by  the  Weather  Facility  for 
more  than  48  hours  In  advance.  Accordingly,  users  should 
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avoid  requesting  plots  for  time  periods  beyond  36  hours  from 
the  time  the  request  Is  submitted.  If  a  plot  Is  requested 
which  is  later  in  time  than  the  last  computed  weather/  the 
information  will  be  provided  for  that  plot  and  a  caution 
statement  will  appear  in  the  message  reply. 

c.  Numbers  of  Initial  and  Datum  Positions 

No  more  than  nine  initial  positions  may  be 
entered/  and  no  more  than  nine  datum  times  will  be 
calculated/  in  any  one  run  of  the  program. 

The  user  Is  encouraged  to  request  whatever 
number  of  datum  points  may  be  useful  to  him.  Thus  if  a 
search  is  planned  for  0800Z/  the  search  planner  may  wish  to 
request  additional  data  for  target  position  at  1000Z/  1200Z/ 
1400Z/  and  1600Z. 

2.  Leeway 

it  is  suggested  that  further  experimentation  be  done 
in  the  area  of  survivor  craft  leeway  to  determine  the 
general  shape  of  the  leeway  function,  it  is  expected  that 
this  shape  will  be  close  to  that  of  the  llferaft  drift 
tables/  and  to  the  formula  of  Witting  given  in  Ref. 3. 

\j  -  k  (35) 

and  which  is  used  in  the  time-sharing  version  of  the 
program  for  all  leeway.  This  same  function  is  used  in  both 
versions  to  approximate  llferaft  leeway. 
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When  better  date  Is  available,  the  program  at  Fleet 
Numerical  can  be  changed  to  accomodate  the  curved  leeway 
function. 

3.  Testing 

It  is  suggested  that  a  rigorous  program  of  testing 
and  evaluation  of  the  results  of  the  Search  Planning  Program 
be  initiated.  If  the  program  proves  to  be  more  accurate 
than  other  methods,  its  contribution  to  the  search  planning 
process,  and  to  oceanographic  research  In  general,  will  be 
substantial . 

4.  Operational  Use  In  the  Coast  Guard 

Although  the  program  will  be  available  to  the  Coast 
Guard  via  the  facilities  of  the  Fleet  Numerical  Weather 
Central,  there  are  further  possible  advantages  of  such  a 
program  which  are  not  available  through  these  facilities. 

a.  Optimum  Search  Areas 

The  computer  method  of  solving  for  optimum 
searches  shown  In  part  IV  of  this  thesis  could  be  added  to 
the  Search  Planning  Program. 

b.  Other  Applications 

Computation  of  on-scene  aircraft  endurance, 
time  enroute  search  areas,  division  of  search  area  among 
different  search  units,  maximum  rescue  coverase  intercepts, 
fixes  from  radio  direction  finder  bearings,  and  many  other 
calculations  familiar  to  search  planners  could  all  be 
available  instantly  from  a  remote  computer  terminal. 
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The  Coast  Guard  Amver  computer  center  at  Governor's 
Island  New  York  regularly  receives  weather  data  by  wire  from 
Fleet  Numerical  Weather  Central  which  is  used  to  produce 
weather  maps.  The  fields  required  for  the  5earch  Planning 
Program  could  also  be  transmitted  to  a  central  Coast  Guard 
computer  facility.  Remote  terminals  connected  to  this 
facility  and  Installed  in  Rescue  Coordination  Centers 
(RCC's)  could  provide  the  output  from  this  program,  Amver 
surface  picture  plots,  and  other  useful  data  directly  to 
search  planners. 

3.  DETERMINATION  OF  SEARCH  AREAS 

It  Is  recommended  that  the  findings  of  parts  IV  and  V  of 
this  thesis  be  evaluated  by  the  Coast  Guard  and  applicable 
portions  be  included  in  the  National  Search  and  Rescue 
Manual . 

C.  FUTURE  RESEARCH 

There  are  many  unexplored  areas  In  the  field  of  Search 
Planning  In  which  further  research  may  yield  valuable 
resul ts. 

1.  The  Trackl Ine  Search 

Although  the  type  of  point-datum  search  presented  In 
parts  IV  and  V  are  computed  In  two  dimensions,  there  are 
times  when  the  target's  distribution  Is  centered  on  a  line 
of  position  (or  trackllne). 


which  presents  a  problem 


essentially  one  dimensional  In  nature. 


If  a  boat  departs  point  A  Intending  to  follow  a 
certain  track  AB  to  point  B,  and  does  not  arrive  at  point  B, 
and  If  no  other  Information  Is  available  as  to  his  probable 
position,  a  one-dlmenslonal  or  trackline  search  problem  is 
presented. 

a.  Partial  Solution 

Assuming  that  It  is  decided  to  begin  the 
searching  efforts  by  executing  just  one  sweep  along  the 
length  of  the  trackline,  what  Is  the  probability  that  the 
target  will  be  located?  The  following  analysis  Is  for  the 
inverse  cube  law  of  detection. 

Given  the  standard  deviation  of  the  target's 
error  of  position  CT,  and  given  sweep  width,  W,  where 


VV  -  2  yj  z  rrrr) 


(36) 


Solving  for  m,  which 
other  factors, 

w' 


m  = 


8  TT 


Includes  visibility  and 


(37) 
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which  has  the  solution 


W 

\Tz  Cr 


(45) 


b.  General  Solution 

The  solution  for  probability  P  of  detecting  a 
target  on  a  trackline  search  for  more  than  one  sweep 
Involves  a  more  complicated  express  Ion,  one  which  does  not 
lend  itself  to  such  a  simple  solution.  Since  a  trackline 
search  is  executed  one  sweep  at  a  time/  a  general  solution 
for  the  optimum  distance  of  each  sweep  from  the  trackline  is 
indicated/  and  should  be  the  subject  of  further  study. 

2.  Optimum  Search  Speed 

The  Search  and  Rescue  Manual  prescribes  a  "moderate" 
speed  for  searching.  Koopman's  formula  for  sweep  width/  from 
( 9  )  gives 


This  indicates  that  sweep  width  Is  Inversely 
proportional  to  the  square  root  of  the  relative  speed  of  the 
target/  w.  When  search  speed  Is  very  large  in  comparison  to 
target  speed/  w  takes  on  the  value  of  the  searcher's  speed. 
As  w  increases/  however/  the  amount  of  search  path  length  L 

certain  period  of  time  increases 

CO 


that  can  be  covered  In  a 
proportionally: 

L= 
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It  follows  that  the  amount  of  search  effort,  W  x  L, 
that  can  be  effected  by  a  certain  number  of  aircraft  ir.  a 
ftxed  period  of  time  increases  as  the  square  root  of  search 
speed : 


J&L 


-  ^ 


W  L 
A 


*2- 


VNL  .  A .  ^ 


<x  lx) 


4,/ 


1  x  hcu 


(47) 


It  may  be  argued,  then,  that  when  time  Is  a  factor, 
such  as  where  persons  in  distress  may  not  survive  for  long 
periods  of  time,  searching  should  be  accomplished  at  the 
highest  speed  practicable.  In  addition,  since  the  coverage 
of  a  designated  searcr  area  will  take  less  time  at  a  higher 
speed,  lookouts  may  encounter  less  fatigue,  and  the 
effectiveness  of  the  search  may  be  enhanced. 


D.  CONCLUSION 

In  conclusion,  there  Is  a  need  for  more  scientific 
research  In  the  practical  aspects  of  search  planning,  and  a 
corresponding  need  for  the  education  of  search  planners  In 
the  fields  of  probability  and  search  theory. 
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APPEND! X  A 


USER'S  GUIDE  TO 
SEARCH  PLANNING  PkUGRAM 


1 .  flfts&tlaiifla 

The  MONTEREY  SEARCH  PLANNING  PROGRAM  Is  a  computer 
program  whi.h  solves  open  ocean  drift  plotting  problems#  and 
provides  datum#  mlnimax  datum.  Total  Probable  Error  of 
Position,  and  Search  Rad  us  for  both  simple  and  complex 
searches.  The  use  of  Mis  program  Is  available  upon 
request,  on  a  test  basis,  to  ail  Rescue  Coordination  Centers 
through  the  facilities  of  the  0.;».  Navy  Fleet  Numerical 
Weather  Facility,  Monterey,  California. 

2.  input  Paid 

a.  Leeway  Code 

For  this  program,  leeway  is  expressed  as  a  two-digit 
code  as  shown  in  Table  1. 
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00  No  Leeway 

10  Liferaft  with  Drogue 

11  Liferaft  w 1 1 hou t  Drogue 

12  Li feraf t  (Mini  ma  x  ) 

13  -  19  Downwind  drift  (mini max) 

Minimum  =  15  of  wind  velocity 
Maximum  -  Second  digit  as  %  of 
Wind  Veloc i ty. 

20,30,610.  Downwind  drift  at  25,31,  etc. 
of  the  wind  velocity 

21  -  99  Minimax  drift  at  some  fixed  percent 

of  wind  velocity  equal  to  first  digit; 
Deflection  to  the  right  and  left 
of  the  wind  line  equal  to 
10  degrees  times  the  second  digit. 

Table  1.  Leeway  Code 


Example:  43  ■  4%  of  wind  velocity 
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♦.  30  degrees  off  wind  line 

b.  Number  of  Starting  Points 

The  number  of  starting  (or  distress)  positions  must 
be  Indicated.  The  maximum  number  that  can  be  supplied  is 
nine. 

c.  Number  of  Datum  Times 

The  number  of  datum  (or  search)  times  for  which 
plots  are  requested  must  be  supplied.  The  maximum  number  of 
these  is  also  n 1 ne. 

d.  Starting  Point  Data 

Starting  positions  must  be  entered  In  chronological 
order  and  each  point  must  be  fully  Identified  by: 

(1)  Date-time  group  (time  zone  Zu’u) 

(2)  Month 

(3)  Latltuoe  and  Longitude  In  degrees  and  minutes. 
Tenths  of  minutes  must  be  specified.  The  number  of  these 
starting  positions  must  agree  with  paragraph  2. a. 

e.  Datum  Times 

Complete  Date-time  group  (including  month)  for  each 
datum  time  must  be  supplied,  In  chronological  order.  The 
number  of  these  datum  times  must  agree  with  paragraph  2.c. 

f.  Last  Known  Position 

This  Is  the  last  known  position  of  the  search  object 
(see  Para.  611,  National  Search  and  Rescue  Manual)  to  be 
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used  In  computing  DR  position  error.  This  position  must  be 
complete  with  degrees,  minutes,  and  tenths  of  minutes.  If 
DR  position  error  is  not  applicable,  enter  one  of  the 
Initial  positions  from  2.d. 

g.  Error  Factors 

These  factors  must  all  be  specified  to  at  least  one 
decimal  place,  in  the  following  order: 

(1)  Probable  navigation  error  of  search  craft. 

(2)  Initial  probable  error  In  position  of  the 

search  object. 

(3)  Percent  of  DR  distance  to  be  applied  as  DR 

error  (See  Para.  611,  National  Search  and  Reset e  Manual). 
Enter  0.0  If  not  applicable. 

h.  Search  Designators 

These  should  be  one-dlglt  numbers,  from  0  to  5, 

which  describe  which  Safety  Factors  are  to  be  applied  to  the 
Total  Probable  Error  of  position  for  each  datum  time.  One 
digit  must  be  supplied  for  each  datum  time. 

0  No  search  radius  is  to  be  computed  for 

that  datum  time 

1,2, 3, 4, 5  The  Safety  Factor  corresponding  to 
search  1,2, 3, 4,  or  5  Is  to  be  applied. 

3 .  Message  Request  Format 

A  message  request  to  Fleet  Numerical  Weather  Central 
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will  insure  that  the  data  is  entered  correctly  into  the 
computer : 

FM  CEA 

TO  FLENUMWEACEN  MONTEREY 

BT 

UNCLAS 

DISTRESS  HOMEY  65  OVEROUE 

1.  REQUEST  SAR  PLOT. 

2 .  OATA : 

A.  12/2/3// 

B.  03|9||Z  AUG, 3521. IN, 05216. 5W/I3120IZ  AUG, 

36M. IN, 05140. |W// 

C.  I316MZ  AUG/P322MZ  AUG/040200Z  AUG// 

0.  35M. IN,  15160. IW/S. I/ll. |/|.  11/112// 

This  sample  message  requests  a  minimax  plot  of  a 
liferaft  from  two  Initial  positions,  with  data  requested  for 
three  different  times  (para.  2. A.).  The  Initial  positions 
and  date-time  groups  are  given  In  para.  2.B.  (note  that  2 
digits  are  used  to  express  degrees  latitude  and  2  digits  are 
used  for  degrees  longitude). 

Paragraph  2.C.  gives  the  date-time  groups  of  the  datum 
times.  Paragraph  2.0.  gives  the  last  known  position,  error 
factors,  and  search  designators.  for  this  example,  the 
first  search  I s  to  be  computed  for  031600Z  and  a  second 
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search  for 

040200Z. 

No  search  radius  is  to 

be  computed 

for 

032200Z; 

pos i t i ons 

and  position  errors 

will  still 

be 

provided. 

4.  Caution;  Offshore  Regions 

This  program  is  designed  primarily  for  open  ocean 
search;  when  data  is  supplied  for  water  areas  within  100 
miles  from  shore,  answers  should  be  checked  by  hand 
calculations.  Requests  for  SAR  Plots  within  10  miles  from 
shore  should  be  avoided. 

Weather  information  used  in  this  program  Is  applicable 
to  the  Northern  Hemisphere  only. 
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